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2 . Abstract 



We use Dirac's constraint dynamics to obtain a Hamiltonian formulation 
of the relativistic iV-body problem in a separable two-body basis in which 
the particles interact pair- wise through scalar and vector interactions. The 
resultant iV-body Hamiltonian is relativistically covariant. It can be easily 
separated in terms of the center-of-mass and the relative motion of any two- 
body subsystem. It can also be separated into an unperturbed Hamiltonian 
with a residual interaction. In a system of two-body composite particles, the 
solutions of the unperturbed Hamiltonian are relativistic two-body internal 
states, each of which can be obtained by solving a relativistic Schrodinger- 
like equation. The resultant two-body wave functions can be used as basis 
states to evaluate reaction matrix elements in the general iV-body problem. 



We prove a relativistic version of the post-prior equivalence which guarantees 
a unique evaluation of the reaction matrix element, independent of the ways 
of separating the Hamiltonian into unperturbed and residual interactions. 
Since an arbitrary reaction matrix element involves composite particles in 
motion, we show explicitly how such matrix elements can be evaluated in 
terms of the wave functions of the composite particles and the relevant Lorentz 
transformations. 
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I. INTRODUCTION 



The theoretical description of the relativistic iV-body problem is an interesting but per- 
plexing problem in nuclear and particle physics. It involves the possibility of relativistic 
bound states of composite particles on the one hand and the reaction matrix elements 
between these particles at high energies on the other hand. The relativistic iV-body prob- 
lem is complicated by the fact that the relativistic bound state problem is basically non- 
perturbative in nature and cannot be solved by the conventional perturbative quantum field 
theory. One also needs to describe the reaction process relativistically by making use of the 
results from the bound states. 

Much progress has been made in the study of relativistic two-body bound state prob- 
lems |I|-16|. In the 1970s, several authors used Dirac's constraint mechanics |l| to attack the 
relativistic two-body problem at its classical roots @], successfully evading the so-call "no 
interaction theorem" |J. The quantum version of the constraint approach was extended to 
pairs of spin one-half particles. The results were two-body quantum bound state equations 
that correct the defects in the Breit equation, correct the defects in the ladder approxima- 
tion to the Bethe-Salpeter equation, and control covariantly the relative time and energy 
variables Those bound state equations for fermions are the two-body Dirac equations 
of constraint dynamics, which we shall also call the constraint equations ||. They possess 
a number of important desirable features, some of which are unique. For example, they re- 
move the deficiencies of earlier approaches in which the spin dependence of the potentials is a 
patchwork of semirelativistic corrections determined by field theory ||. In contrast, the spin 
dependence in constraint dynamics is determined naturally by the Dirac equation structure 
of each of the constraint equations together with the assumed covariant interactions |j[[7|. 
The equations of constraint dynamics are manifestly covariant while yielding simple three- 
dimensional Schrodinger-type equations, like those of their nonrelativistic counterparts ||[7| . 
(This particular feature will ultimately be of crucial importance in the discussions given in 
this paper.) These constraint equations have passed numerous tests showing that they repro- 
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duce correct QED perturbative results when solved nonperturbatively 0. In addition, the 
Dirac forms of these equations automatically make it unnecessary to introduce ad hoc cutoff 
parameters, which are needed in most other approaches to regularize singular potentials. 
The relativistic potentials appearing in the constraint equations are related directly to the 
interactions of perturbative quantum field theory. In non-perturbative QCD as applied to 
meson spectroscopy, they may be introduced phenomenologically |10| and can be regarded 
as an anticipation of potentials that may eventually emerge from lattice gauge theory. 

In the present paper we will extend the constraint equations beyond the scope of the 
two-body problem and will, for simplicity, limit ourselves to the constraint description of 
spinless particles. We first summarize the constraint approach for two spinless particles. For 
each particle, one ascribes a generalized mass shell constraint which includes the interaction. 
The constraints must be consistent with each other and this in turn restricts the dependence 
of the interactions on the relative coordinates, eliminating both the relative time and energy 
variables in the CM system. The resultant equations correspond to a Bethe-Salpeter equa- 
tion whose kernel and Green's functions are constrained by the requirement of P-q =0, where 
P and q are the total and relative momenta of the two-body system ||. In particular, from 
the Bethe-Salpeter equation with this constraint, one can derive pH|J7[] the "quasipotential 
equation" of Todorov |§ , which is a Schrodinger-like inhomogeneous integral equation where 
the quasipotential $ is related to the scattering amplitude in perturbative quantum field 
theory. Other methods of reducing the Bethe-Salpeter equation have also been suggested 

mm 

At present, the relativistic treatment of the iV-body problem and the reaction of com- 
posite particles at high energies have not advanced as much as in the relativistic treatment 
of two-body bound states. An investigation of the relativistic iV-body problem was previ- 



ously carried out by Sazdjian flllf . In contrast, the non-relativistic iV-body problem and the 
non-relativistic description of the reaction of composite particles have been well developed 
|T7^pU||. One has, for example, the distorted- wave Born approximation method in nuclear 
reactions jjlTj] and the quark-interchange model in hadron reactions [|l8l for reactions between 



composite particles. One first solves for the wave functions and determines the interaction 
between the constituents using the energy levels of bound states. Then one uses the same set 
of interactions and wave functions to calculate the reaction matrix elements for the reaction 
of composite particles. There is also the potentially different values of the post or prior forms 
of the reaction matrix element, which distinguish whether the interaction occurs before the 
rearrangement or after the rearrangement of the constituents. The post-prior equality of 
the reaction matrix elements is attained when the interaction and wave functions that are 
used in the calculation of the overlap integral are the same as the interaction and the wave 
functions obtained in the bound state analysis. 

The formalism we shall develop can be applied to many processes of interest. For exam- 
ple, in high-energy heavy-ion collisions, the investigation of the dynamics and the properties 
of the produced hadron matter involve the reaction cross sections and the reaction matrix 
elements between the produced hadrons at relativistic energies. Most of these cross sections 
and reaction matrix elements cannot be measured experimentally. A reliable theoretical 
relativistic model which describes the iV-body problem of the constituents is needed for 
their evaluation. It is necessary to generalize the nonrelativistic reaction model to study 
reactions between composite particles at relativistic energies. It is also of interest to study 
the relativistic iV-body problem in order to investigate an assembly of composite particles 
and their clustering or molecular states. 

We shall take advantage of previous advances in the understanding of the relativistic 
two-body problem in our present study of the relativistic iV-body problem. In this work, we 
specialize in iV-body dynamics with just pair-wise interactions between particles. We shall 
write down the Hamiltonian formulation of the relativistic iV-body problem which allows 
an easy separation in terms of the center-of-mass and the relative motion for any two-body 
system and a simple separation of the unperturbed Hamiltonian and the residual interaction. 
We can then use the relativistic two-body bound states as basis states for the investigation of 
the relativistic iV -body problem. We can construct a proof of the "post-prior" equivalence 
in relativistic dynamics, which guarantees that the reaction matrix element is independent of 
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the different ways of partitioning the unperturbed Hamiltonian and the residual interaction 
from the iV— body Hamiltonian. 

This paper is organized as follows. In Section II we review important aspects of the 
constraint approach and introduce the relativistic two-body Hamiltonian and bound state 
eigenvalue equation. We show that the relativistic two-body solution in the CM rest system 
has the simplicity of its nonrelativistic counterparts. In Section III we generalize the two- 
body Hamiltonian to the case of N particles. We discuss some applications of the present 
formulation in Section IV. We show how to utilize the constituent two-body wave functions 
for the iV-body problem. Because of the analogy to their nonrelativistic counterparts, these 
states can be used as basis states to evaluate a general reaction matrix element in the 
general iV-body problem. We prove the relativistic version of the post-prior equivalence for 
the reaction matrix elements. As an illustration and a problem of practical interest, we 
consider in Section V the reaction of four particles in two composite systems where the non- 



relativistic treatment has already been formulated by Barnes and Swanson [ I8| . In Section 
VI, we explicitly obtain the wave function in second-quantized form so as to construct the 
overlap integral for the reaction matrix element. Such a reaction matrix element involves 
states of composite particles in motion. Thus we show in Section VII explicitly how to 
evaluate such elements in terms of the wave functions of each composite particle, and develop 
the relevant Lorentz transformation laws required. Section VIII summarizes and points to 
future problems. 

II. HAMILTONIAN FORMULATION OF THE 2-BODY PROBLEM FROM 

CONSTRAINT DYNAMICS 

We can formulate the relativistic treatment of the two-body problem for spinless particles 
pT| , |22|j in a way that has the simplicity of the ordinary non-relativistic two-body Schrodinger 



equation and yet maintains relativistic covariance. Including spin and generalizing to differ- 
ent types of interactions can be carried out in a more complete framework @,|J. 



For each particle we assume a generalized mass shell constraint of the form 

Wi|V>) = for £ = 1, 2 (2.1) 

where 

Hi^pl-ml-^, (2.2) 

and $i and $2 are two-body interactions dependent on x%2- One constructs the total Hamil- 
tonian 7i from these constraints by 

H = XiHi + A 2 W 2 , (2.3) 

(with A, as Lagrange multipliers). In order that each of these constraints be conserved in 
time we must have 

{H l ,HM = i^\^) = 0. (2.4) 

As a consequence, we have 



[Hi, \iH\ + 

{[^,A 1 ]W 1 |^)+A 1 [W i ,W 1 ]|^) + [W i ,A 2 ]W 2 |V') + A 2 [W i ,W 2 ]}|^)=a (2.5) 

Using Eq. (|2j]) , the above equation leads to the compatibility condition between the two 
constraints, 

[HuKM =0. (2.6) 
Since the mass commutes with the operators, this implies 

([p?,$ 2 ] + [$ 1 ,^] + [$ 1 ,$ 2 ])|^) = 0. (2.7) 
The simplest way to satisfy the above equation is to take 

$1 = $2 = (2.8) 



which is a kind of relativistic Newton's third law. Here, the transverse coordinate is defined 
by 

x u ± = x^irj^ - P^Pu/P 2 ), (2.9) 
where P is the total momentum 

P = Pi+P2- (2.10) 
The choice of the two-body potential Eq. ( |2.8|) leads to 

[Hi,n 2 ]\i>) = ip ■ d xia $(x ± )\if>) = o, (2.ii) 



and the compatibility condition ( |2.6| ) is satisfied. 

The two-body Hamiltonian 7i determines the dynamics of the two-body system. Its 
equation of motion is 

W|V>) = 0. (2.12) 

This equation describes both the center-of-mass motion and the internal relative motion. 
To characterize the center-of-mass motion, we note that since the potential $ depends only 
on the difference of the two coordinates we have 

[P,HM = 0. (2.13) 

(This does not require that [P, Aj] = since the TCi\ip) = 0.) Thus, P is a constant of motion 
and we can take to be an eigenstate state characterized by a total momentum P. 

To separate out the internal relative motion from the center-of-mass motion, we introduce 
the relative momentum q defined by 

Pi • P 

Pi = ^prP + q, (2.14) 

P2 = ^fP-q, (2.15) 



where the first term on the right hand side of the above two equations is the projection of 
each momentum onto the total momentum. The above definition of the relative momentum 
guarantees the orthogonality of the total momentum and the relative momentum, 

P-q = 0, (2.16) 



which follows from taking the scalar product of either equation with P. From Eqs. Q2.14 ) 
and ( 2.15| ) this relative momentum can be written in terms of p\ and p 2 as 



I p2 L \JP 2 

where 

pi ■ p _ p 2 + p\ - p\ 



1 Pl -^L= P2 (2.17) 



P 2 2VP 2 

e, = v -^L = ^±4zA (2 . 18) 
VP 2 2VP 2 



are the projections of the momenta p\ and P2 along the direction of the total momentum P. 
Using Eqs. ( [2.1| ) and ( |2.8| ) and taking the difference of the two constraints, we obtain 

{pl-plm = {m\-ml)\i,). (2.19) 

Thus on these states \ip) we have 

P 2 + m\ — m\ 

E\ = ■ 



2VP 2 

P 2 + ml — m? , 

£2 = A- -■ 2-20 

2^/P 2 



Using Eqs. ( |2.14j ), ( |2.15| ), and Eq. ( |2.16j ), we can write 7i in terms of P and q: 

U\i> ) = {Ai[e? - m\ + q 2 - $(z ± )] + \ 2 {el -m 2 2 + q 2 - $(x±)]}|V> 

= (Ax + A 2 )[6 2 (P 2 ; m\,m\) + q 2 - $(x ± )]\ip) = 0, (2.21) 

where 

h 2 {P 2 ,m\,m\)=e\-m\=el-m\ = ^(P 4 - 2P 2 (m\ + m\) + [m\ - m\f). (2.22) 



Equation ( [2.2 1|) contains both the center-of-mass momentum P [through the 
b 2 (P 2 , m\, wij)] and the relative momentum q. This constraint equation of P and q can 
then be solved by the method of the separation of variables. That is, we introduce the 
bound state eigenvalue M to separate Eq. ( p.21| ) into the following two equations for the 



center-of-mass motion and the internal motion 

{p 2 -M 2 }|^> = 0, (2.23) 

and 

(Ai + A 2 ) [q 2 - + b 2 (M 2 , m\, m 2 2 )} \f) = 0, (2.24) 

where we have used the first equation on the eigenstate so that b 2 {P 2 ,m\,m\) becomes 
the standard triangle function indicative of the presence of exact relativistic two-body kine- 
matics: 

b 2 {M 2 , ml, m 2 2 ) = Jj^ {M 4 - 1M 2 {m\ + ™l) + {m\ - m 2 2 ) 2 } . (2.25) 

The eigenvalue equation Eq. ( |2.24| ) for the relative motion is independent of the Lagrange 
multipliers. It is nonetheless convenient to choose Aj = l/2mj so that the resultant 
Schrodinger equation matches the non-relativistic two-body Schrodinger equation term by 
term. Such a choice also helps us obtain useful simplifications in the relativistic iV-body 
problem in later sections. In particular, the iV-body Hamiltonian can be easily separated 
into pairs of two-body Hamiltonians. This separation makes it easy to introduce the unper- 
turbed Hamiltonian and residual interactions. 

We note that because of the orthogonality of P and q, we can write Eq. ( [2.24j ) in the 
form 

(i + {q2± ~ $(;E±) + b2 ( M ^ m l m D}\^) = °> ( 2 - 26 ) 

where q± = q — q ■ PPj P 2 = q. [Note that if the relative momentum were defined in terms 
of Eq. (p.20|) instead of Eq. ( p,18|) then we would have 
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q-P\ip) = (2.27) 

but not q ■ P = so that q 2 \ip) — q\\ip)- In either case the coefficients £j are invariant and 
hence Eq. ( |2.17| ) the same form regardless of which frame it is evaluated in.] 

We show below how the eigenvalue M is related to the eigenvalue obtained in a non- 
relativistic Schrodinger equation. We go to the center-of-momentum system where q = 
q±(0,q) and x± = (0, r) (relative energy and time thus being removed from the problem). 
We then have the equation for the relative motion, 

f q 2 $(r) b 2 ) , v . . 

{l + ^-^j 1 ^ 0, (2 - 28) 

where \i is the non-relativistic reduced mass, 

a = 2.29 

mi + mi 

We can cast Eq. ( |2.28|) into the usual form of a non-relativistic Schrodinger equation. By 



renaming $/2/z as V12, and 6 2 /2/i as E, Eq. ( p.24|) becomes 



f^+^|V> = ^>. (2-30) 

The above Schrodinger equation can be solved to give the eigenvalue E. Then, from the 
equation 6 2 (M 2 ,m 2 ,m2) = 2fiE, one can solve for M in terms of E and obtain 



M = y 2fiE + m 2 + y 2fiE + m\. (2.31) 

It is easy to show from this that in the limit of very weak binding, the nonrelativistic limit, 
we have the familiar result 

M = m 1 + m 2 + E. (2.32) 

If one is only interested in the effect of exact two-body relativistic kinematics with V12 
an energy-independent nonrelativistic potential, the bound state eigenvalue M for the rel- 
ativistic two-body problem is related to the eigenvalue E of the nonrelativistic problem by 
Eq. ( [2.31|) . It is important to note, however, that the potential Vu in relativistic constraint 
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dynamics includes relativistic dynamical corrections as well. These corrections include de- 
pendences of the potential on the CM energy M and on the nature of the interaction. For 



spinless particles interacting by way of a world scalar interaction S, one finds f22],|21 

_$ = 2m M S 
2/i ~~ 2/i 



V a = * = 2m " S + S2 (2.33) 



where 



M 

while for (time-like) vector interaction A, one finds P J22| , |24}| 



mu = - J r 7 1 , (2.34) 



2/i 2/x V ' 

where 

M 2 — m? — ml , 

£ " = 2jt~^ (2 - 36) 

and for combined space-like and time-like vector interactions (that reproduce the correct 
energy spectrum for scalar QED [|2"Tf ) 

$ _ 2e M A - A 2 + V 2 log(l - 2A/M) 1 ' 2 + [Vlog(l - 2A/M) 1 / 2 ] 2 
12 2/i 2/x ' 

The variables mu and £m (which both approach \i in the nonrelativistic limit) were intro- 
duced by Todorov in his quasipotential approach and are called the relativistic reduced 
mass and energy of the fictitious particle of relative motion. In the nonrelativistic limit, 
$ approaches 2/j,(S + A). In the relativistic case, the dynamical corrections to Vu referred 
to above include both quadratic additions to S and A as well as CM energy dependence 
through rriM and Em- This latter point implies that the effective potential Vi 2 depends on 
the eigenvalue E (or M) to be evaluated. One can obtain the mass of the bound state M 
by an iterative procedure. One starts with an estimated M (or E) value and obtains the 
potential Vu- Equations ( |2.3(J| ) and ( |2.31| ) can then be used iteratively to obtain successively 



improved values of V12 and the eigenvalue M (or E). 
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We note in passing that since 

b 2 = e 2 M - m 2 M , (2.38) 

we can write the Schrodinger-like equation for combined scalar and (time-like) vector inter- 
actions as 

{q 2 + (m M + S) 2 - {e M - A) 2 } = 0, (2.39) 

which is suggestive of a Klein-Gordon equation for an effective particle of relative motion. 
This bound state equation incorporates not only the correct relativistic kinematics but also 
the correct relativistic dynamical corrections through order 1/c 2 and higher, depending 
on the input. It does it without the necessity of introducing complicated momentum- 
dependent Darwin-like interactions, thereby retaining the simplicity of the non-relativistic 
Schrodinger equation. Furthermore, the potentials in these equations are connected to those 
of Wheeler- Feynman electrodynamics (and its scalar counterpart) [[2^,[2(J. They have been 
obtained systematically from perturbative quantum field theory |7],[],|II]] and from an eikonal 
summation of Feynman diagrams [[27] . 

In summary, Eqs. ( |2.28[ ), ( |2.30| ), and ( |2.31| ) provide a useful way to obtain the solution of 
the relativistic two-body problem for spinless particles in scalar and vector interactions. In 
other works they have been extended to include spin and have been found to give an excellent 
account of the bound state spectrum of both light and heavy mesons using reasonable input 
quark potentials (Tyj2^ . 



III. HAMILTON! AN FORMULATION OF THE iV-BODY PROBLEM FROM 
CONSTRAINT DYNAMICS: SEPARABLE TWO-BODY BASIS 

The above treatment of the two-body problem can, to some extent, be generalized to the 



case of the iV-body problem [IT]. Our approach differs and extends the work of Sazdjian 



in that we choose to formulate the iV-body problem in a separable two-body basis. 
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We consider a system of TV particles. For each particle, we specify a generalized mass 
shell constraint of the form 



7U\i>) = for i = l,..,N (3.1) 

where 

N 

W« = p?-m?- X) *a-Wi. (3.2) 

The $jj are two-body interactions dependent on a;^, and PVj are possible iV-body forces 
(JV > 2). We construct the total Hamiltonian for the system as 

N 

H=J2 X i n i ( 3 - 3 ) 

i=i 

where A^s are the Lagrange multipliers. For each constraint to be conserved in time, we 
must have 

{H^HW = li-^HtW = 0. (3.4) 
From Eq. ( |3.3|) and Eq. ( |3.1|) , we must have the compatibility condition 

[HitHjM = 0. (3.5) 
Now we attempt to expand out the above equation. For a fixed pair of % and j, we have 

N N N N 

(-[pi E - [ E <w?] + [E + Wi, E + 

k,kj^j k,ky^i k,k^j 

-\piWi\-[w hP )\m 

= 0. (3.6) 
Motivated by the form of our two-body solution, we assume that 

$jj $ ' (3.7) 

in which 

{%ij±)v — x ij V^u — (Pij) n{P%j)v I 'Pfj i (3-8) 
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and Pij = pi + pj. This implies that the iV-body forces must be present and satisfy 

N N 

([£ + W u £ ^ + Wj] - \pl W 3 ] - [W h p)\m = 0. (3.9) 

I Ij^i k k^j 

These are very complicated equations, and, unlike the two-body case, there are no known 
closed- formed solutions [p| , ^| . Evidently these forces are dependent on the two-body forces 



themselves |p9|j . In practice, one often ignores these many-body forces and considers only 
pair-wise interactions, as we will do in our subsequent computations of reaction matrix 
elements in Section VII. That is, under the approximation in which we set Wi = 0, we 
can view the particles as interacting with each other via two-body interactions in a pair- 
wise manner. However, in most of our formal analysis in this paper, we will retain these 
many-body interactions. 

The conservation of the constraints in time depends only upon the compatibility of 
these constraints and does not depend on the choice of the Lagrange multipliers. This 
arbitrariness in the choice of the A, is similar to a kind of gauge invariance. Choosing 
a particular set of A, is analogous to choosing a gauge and can be done for convenience. 
In the general iV— body formalism we will find it convenient to choose Aj = l/2mj. Such 
a choice has many advantages. First it leads to a simple correspondence with the non- 
relativistic two-body and iV-body Hamiltonians. Secondly this choice depends only on the 
particle in question and not what other particle it is linked with. In the general iV-body 
formalism we will find it convenient not to have a preferred pairing of 2-body composite 
subsystems. The choice Aj = mj/2 avoids this. Finally, this allows the TV-body Hamiltonian 
to be conveniently separated into a nonperturbative part and the residual interactions part. 
Such a correspondence helps one generalize the post-prior equivalence of reaction matrix 
elements from the non-relativistic case to the relativistic case. 

With this choice, the relativistic iV-body Hamiltonian is 

N 1 N N <f>.. w 

i=l i=i j j>i Z PU i 

N I N N W 
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where fly = rnirrij/ (mi+rrij), and we have introduced the simplified notation = $jj/2/ijj. 
(Note that in light of the above forms of Eqs. ( |2.33| ) and ( p.36| ) for <J>j, the choice A, = l/2mj 



gives the correct nonrelativistic limit for V%j.) The dynamics of the relativistic A^-body 
system is determined by the search for the state such that 

i=i 1 i=i j>i i 1 

Even without the many-body forces Wi , this equation is very difficult to solve because the 
potentials depend on the momenta = Pi+Pj (through x^jj which are not constants of 
motion for the N— body system; forces outside of the ij system produce time-dependent P^-. 
On the other hand, if one uses the two-body Hamiltonians to generate basis states, then for 
those states and Hamiltonians, one can regard the P^ as constants of motion. That is, the 
Hamiltonian can be separated in terms of two-body Hamiltonians plus residual interactions 
regarded as perturbations. This greatly simplifies the problem. 

IV. SOME APPLICATIONS 



The quadratic form of the momentum operators pi in the iV-body Hamiltonian Eq. (|3. 1 1 



makes it easy to manipulate the momentum terms to obtain the center-of-mass momentum 
and other relative momenta. The potential term in the equation appears in a way similar 
to that in which it appears in the non-relativistic case. 

The relativistic A^-body equation can be compared to the non-relativistic A^-body equa- 
tion. Introducing e» = p\ — m,, we have 

' (p?-rn?) = «-4^- ^ 



2rrii 1 ' 2rrii 
Equation (|3.11| ) becomes 

W) = {pnr - E - E E v* - E f- 1 } Wi = o, (4.2) 

where Pnr = J2i e i- This is identical to the non-relativistic A^-body Hamiltonian with 
eigenvalue Pnr when \e\i « m» and the Wi are neglected. 
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In the next example, we can examine a system of an even number of iV particles forming 
N/2 composite particles, as in a system of N/2 mesons. For such a system, one can consider 
an initial state of the form 



l^o) = (i2, h), -(ijv/2j'jv/2)}), (4-3) 

in which particles i a and j a form a composite two-body subsystem (i a j a )- Subsequent 
dynamics of the system is determined by the evolution operator containing Ti and the 
reaction matrix element {ip' a \H\ip a ) where 

I^) = |{(<jO,(4,j 2 ),---(4/ 2 jV)})- ( 4 - 4 ) 

For the evaluation of the element (ip' a \TC\ip a ) of the Hamiltonian matrix, the iV-body Hamil- 
tonian can be separated into an unperturbed Hamiltonian Ho and a residual interaction 
Vi, 

H = H + V t , (4.5) 
where the unperturbed Hamiltonian 7i is 

7~Lo = Hhn + Wi 2 32 + •• + ^jv/ 2 ijv/2' ( 4 -6) 

with 

and Vi, the 'prior' form of the residual interaction, is 

N N ' W 

^ = -EE^-E^ (4-8) 

i=lj,j>i i zm i 

where the summation X)£Li Sjj>/ * s carried out with the set {ij} different from those of the 
composite particles {hji},{i 2 j2}, {hjs}, -, and {in^Jn/t}- 



N N ' N N 

EE =EE 

* =1 i.i> 1 »=i j,j>i 



(4.9) 

{lj}^{llj'l},{l2J2},-{»JV/2jjV/2} 
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Then, since 

Wo|{(iljl), (*2, j 2 ), -(^/2jiV/2)}) = 0, (4.10) 

the transition matrix element (■^I'^lV'a) becomes 

mnm = mvi^a). (4.11) 

The separation of W into 7i and V/ is not unique, and it is important to show that 
the reaction matrix element of the residual interaction is independent of the different 
ways of separating out the unperturbed Hamiltonian and the residual interaction. We 
can alternatively choose the unperturbed Hamiltonian to be associated with the state 

Wa) = {{(AjtiA^fi),- -(ijv/^/a)}) such that 

Kimti), -(4/2^/2)}) = o. ( 4 - 12 ) 

associating with the "post" form of the residual interaction V/, 

n = n' + v;. (4.13) 

The quantities 7i' and VJ are defined in a way similar to that given in Eqs. ( |4.6H4.9| ). The 
reaction matrix element between the basis states is 

= mV^a), (4-14) 

which indicates that the reaction matrix element is the same for the "prior" form or the 
"post" form of the residual interactions. It is independent of the way in which we split 
up the total Hamiltonian. This "post-prior" equivalence guarantees the uniqueness of the 
reaction matrix element and insures the usefulness of the perturbation expansion. 

Using this method of separating the total Hamiltonian, the matrix element of 7i between 
any two basis states can be evaluated. The construction of the matrix of the total Hamil- 
tonian will allow one to construct the evolution operator and to follow the dynamics of the 
system. 
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V. SCATTERING OF TWO COMPOSITE PARTICLES 



As another explicit example, we apply our formalism to a problem of practical interest. 
Consider four particles with masses m; and momentum Pi where i — 1,2, 3, 4. We have the 
total Hamiltonian 

«=e^-"B-i:ev^ p.i) 

i=i * i=i j,j>i i=i < 

Using the method of the first section, we can solve for the bound states of mass My for the 
motion of particles i and j interacting with the interaction V%j, 



1 , , 1 



ij 



m = 0. (5.2) 



3 

We can consider the reaction of two composite particles A{12) and 5(34) where particles 
1 and 3 are particles and particles 2 and 4 are antiparticles as in meson-meson scattering, 
interacting through a pair- wise interaction Vij. We study the relativistic quark-interchange 
reaction 

A(12) + 5(34) -> C(14) + £>(32), (5.3) 



as a generalization of the non- relativistic case investigated by Barnes and Swanson |L8| . For 
this reaction with a momentum transfer t 

t = (A - Cf = m\ + m 2 c - 2A C + 2A • C, (5.4) 



the differential cross section in the first-Born approximation is given by [18 



da 1 % 



2 



\M fl \ 2 (5.5) 



dt 64vrs |p^| 2 

where M.fi is 



M fi = {2irypE c A 2E c B 2E c c 2E c D h fi , (5.6) 

and hfi is the reaction matrix element from the initial state A(12)S(34) to the final state 
C(14)D(32) initiated by the residual interaction Vj. In the above equation, a kinematic vari- 
able with the superscript c refers to that variable evaluated in the center-of-mass (collider) 
system. From the above result, hfi has the dimension of l/(mass) 2 . 
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To obtain the reaction matrix element hfi in our relativistic formulation, we need to 
split the total Hamiltonian into the nonperturbative part 7i and the residual interaction 
Vj. This can be carried out in two different ways. In the "prior" form, it is split as 

4 TT7 * 

n = n 12 + n 3i - v 13 - v u - v 23 - vk - E -A (5.7) 

where the unperturbed Hamiltonian is 

^o(prior) = H 12 + H U , (5.8) 

and the residual interaction is 

4 W 

VKprior) = ~V 13 - V 14 - V 23 - V 24 - 7T~- ( 5 - 9 ) 

i=l Zm i 



The reaction matrix element is 



4 



2n5\P A + P b -P c - P D )h fi (pnor) = -(^14^23^13 + V u + V 23 + V 24 + Y / 

i=l 

(5.10) 

In graphic form, if we represent the interaction Vij by a curly line, the first four terms in the 
above matrix element are represented by the four diagrams in Fig. 1. The interaction takes 
place before the rearrangement of the constituents. 



a . . . c 

A ^fa " " '« C 



B 




D 




CI C2 





T2 



Fig. 1. 'Prior' diagrams for the reaction A+B — > C+D. 

On the other hand, if we use the "post" form of splitting the total Hamiltonian, we have 

4 W 

h = h 14 + n 32 - v 13 - v 12 - v 43 - y 42 - E 7T^- ( 5 - n ) 



1 2mj 
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The unperturbed Hamiltonian is 

n (post) = n l4 + n 32 , (5.12) 

and the residual interaction is 

Vr(post) = -V 13 - V 12 - V 43 - V 42 - E T^- ( 5 - 13 ) 

i=i * 

The reaction matrix element is 

4 W- 

2Tl5\P A + Pb-P C - P D )h fi (pOSt) = -(^14^23^12 + Vxa + V42 + V43 + £ 7^1^12^34). 

i=l 

(5.14) 

In graphic form, the first four terms in the above matrix element are represented by the four 
diagrams in Fig. 2. The interaction takes place after the rearrangement of the constituents. 




Tl T2 

Fig. 2. 'Post' diagrams for the reaction A+B — > C+D . 

Therefore, if we start with the prior expression for the matrix element, we have Hvi^vi = 
and 7^34^34 = 0, and we have (cancelling out the iV > 2-body potentials) 

(^14^231^13 + Vu + V23 + ^241^12^34) = (V'mV^I - ^12 - 7^34 + V 13 + V U + V 23 + V 2A 1^12^34) 

= (V'uV'aal - ft 14 - ^23 + Via + + V 42 + V 43 \1p12ip34), (5-15) 

where we have used Eq. fl5.1| ) to write out the Hamiltonian for the two-body system. Because 
TCuI^ia) = and 7^23 1^23) = 0, we have then 

1 

(^14^231^3 + Vl4 + V23 + V 24 \ijj 12 ^ 3i ) = (-014-023 1 V12 + Via + V* 42 + ^^12^34), (5.16) 
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which leads to the relativistic generalization of the post-prior equivalence of the reaction 
matrix element hfi, 

hfi (prior) = h fi (post). (5.17) 

Just as in non-relativistic reaction theory P0 |, the equivalence is possible only when one 
uses the same internal relative wave function for the composite particles in their scattering 
process as in the bound state problem for the individual composite particles. The equivalence 
allows a unique determination of the reaction cross section in the first-Born approximation. 



VI. SECOND QUANTIZATION OF PARTICLES IN A BOUND STATE 

A practical problem arises when one attempts to evaluate the reaction matrix element 

(V>.IW«>, 

N N ' rrr. 

(A\ViW a ) = ((I,*), (idz), ...| -EE - E^K^'D, ■••>, (6-i) 

i=i jj>i i zm i 

where {(ziji), (^2^2), ••■} represents composite two-body subsystems. For a pair-wise inter- 
action Vij in the above equation, we have 

<^a| W a ) = (W, UP'WMia), UP)), (6-2) 

where {a'/3'} is a permutation of {a/3}. The composite wave functions ipi a are usually 
computed in the CM of the (ia) composite particle system (the qiq a system in our exam- 
ple of a system of mesons). However, this is not the same as the so-called collider frame 
(the CM frame of the (ia)-(jf3) meson-meson system) either in the bra or ket states. For 
non-relativistic reactions, the relative wave function of a composite system in the collider 
frame is obtained from the wave function for the composite particle at rest by a Galilean 
boost, and they are related by a simple momentum shift. In the relativistic case, a Lorentz 
boost is needed in place of a Galilean boost. We need, therefore, to discuss the Lorentz 
transformation of the state of the composite system. 

We represent the state of the composite two-body system with a 4-momentum P by 
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(P°|(12) P > = Ae^ J d 4 Pl d A p 2 5\ Pl +p 2 - P)0(p°)0(p°) 

x5(p 2 -m\- $(x ± ))6(pl -m\- $(x ± ))ip(p 1 ,p 2 )\pip 2 ), (6.3) 

where we use the same symbol p\ and p 2 to denote c-numbers and operators, using the 
context to distinguish between them. The factor of 4e 2 e 2 is included so that we obtain 
the usual results in the nonrelativistic limit. The delta function containing the composite 
particle energy P° arises from the projection of the energy eigenstate state \P°) onto the 
state vector |(12) P ) where P = {E,P} and E = VP 2 + M 2 (see also Eq. ( ggg) below). 
Hence, we use the notation (P°|(12)^p) to denote the projection of the state vector |(12)p) 
onto the energy eigenstate state |_P ) • The above state is constructed in analogy to the two 
free-particle states 

(P°|(12) P ) = 4 £l £ 2 J d%d%5 4 ( Pl +p 2 - P)e(p 1 )6(p 2 )5(p 2 1 - m\)5{p 2 2 - mj)^( Pl ,p 2 )\p lP2 }. 

(6.4) 

However, unlike the free-particle state, the state ( |6.3j ) satisfies the simultaneous constraint 
conditions of Eq. ( |2.1| ) 

Hi\(12) P ) = [p 2 -m 2 -^(x ± )]\(12) P ) = for z = l,2, (6.5) 

so that neither particle is on mass shell. The definition ( |6.3| ) also introduces a momen- 
tum space wave function i^( P i,p 2 ) defined so that it has positive constituent energies. We 
emphasize that the momentum eigenstates \ P i P2 ) in this expansion are off shell. That is 

(p 2 i-m 2 )\p lP2 ) t^O. (6.6) 

Our first step is to show that the above bound state composite is a sharp state being zero 
unless P 2 = M 2 where M is the meson bound state mass. By using the total and relative 
momentum operators and Eq. ( |2.2Up (so that Eq. Q2.27D is satisfied and not P ■ q = 0), the 
product of the two delta functions can be written as 

S(pj -m\ - &(x ± ))5(pI -m 2 ,- $(x±)) 
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= 5 (^T + ? - m i " $ (^) + -^-^) 5 (K - m\ y — 2 +2P-q) 
^ ^ +W .^ + i p.J (6.7) 

where 7~£ q = q 2 — +6 2 . We assume that the momentum space wave function ip(pi,P2) 

is an eigenfunction of 7Y g , so that H g i>(pi,P2) = 0. Thus, on such states the above equation 
becomes 

M 2 

5{pi -ml- <S>(x ± ))5(p 2 2 -ml- ${x±)) = 5{P 2 - M 2 )5{P ■ q). (6.8) 

2eiE2 



This shows, as anticipated, that the state defined in Eq. (|6.3|) is sharp and would satisfy 
Eq. (ET23D in addition to Eq. Q. 

We allow the delta function arguments to operate on the momentum states, and use the 



positive energy condition. Then from Eq. (2.17) we have 



0{pl)e{p 2 )8{p\ -m\- Hx ± ))6(p 2 2 -mj- ${x ± )) 
M 2 E E 

= 0(pS)«« ii ^*(p? ~ m £i ~ P ' q)6{P2 ~M £2 + P ' q) ' (6 ' 9) 

where E = VP 2 + M 2 . We would like to express this in terms of a Lorentz transformation 
from the CM system (in which the meson has a mass M). 

When the composite particle is boosted by A to a momentum P, the meson's velocity is 

V = |. (6.10) 

Thus with u = V = P/|P|, the components of A are 

E P'P/% 

Afc = S ik + (7 - l)uiU k = 5 ik + (— - 1) 



M ' P 



2 



Ao = A° = ^v/(7 2 -l) = ^, 

An = 7 = —. (6.11) 

We use the notation p\ = (ex, p*) to represent the four-momenta of the zth constituent in 
the composite particle rest frame, and pi = piA = (siA, P«a) to represent the 4-momentum of 
the ith constituent in the frame boosted by A. Then we have 
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Ef- P • n* 

p» = pi = (Aplf = 1P f + u ■ p* = J± + = e lA , (i = 1, 2) (6.12) 



and 



P E PP • D* 

Pl = p tA = — £i + p* + (-- 1)—^ , (* = 1, 2). (6.13) 



Because p^ + p 2 = we have therefore 



M 

Furthermore, using Eq . (|2 . 1 7| ) we have 



v \ + p ° 2 = VP 2 + M 2 = E, (6.14) 



Pi + P2 = ^1 + ^) = P. (6.15) 



P-q=— Pq* = — P.p? = -— Pp* (6.16) 

M M M y 1 



Hence, Eq. ( |6.9|) can be rewritten as 



9{pl)d{pl)5{p\ -m\- $(x ± ))S(p 2 2 -m\- $(x ± )) 
M 2 

= ~ ( A PI)°)5(P°2 - (A?;) ). (6.17) 

We have therefore 

M 2 r 

(p°i(i2) EP ) = — j d* Pl d 4 p 2 5\p- Pl -p 2 )5(Pi - m)°np°2 - (Aps)>(pi.ft)ipi.ft>- 

(6.18) 

Now the pi and p 2 part of the d A p\ and d A p 2 can be integrated out and the result is 
M 2 r 

(P°\(12) EP ) = —5{P»-p\-p§ j d Pl dp 2 6 3 (P - Pl - p 2 )^(p?Pi,P§P2)|p?Pi,rfp 2 >, 

(6.19) 

where p° = (Ap*)° and p® +p9, = E. In the CM system, p\ + p 2 = M and the state vector is 
(P°\(12) MO ) = d(P -M) J dpldp^ipl + p^Mie^e.pD^ple.p;). (6.20) 
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We introduce the notation |M(P)) denned as 



|M(P)> = J d Pl d P2 5 3 (P - Pl - P2)V(PlPl,P2P2)blPl:P2P2> 



(6.21) 



so that, since in a general frame 



Pl +P2 



= VP 2 + M 2 = E, 



(p°m) EP ) 



M 2 
IP 



5(P { 



- VP 2 + M 2 )|M(P)). 



(6.22) 



The projection of the state vector |(12)^p) onto the time component of the center of mass 



From the above results, the energies of the constituents in a composite particle take 
on fixed values (p®* = si) in the center-of-mass system, while their off-shell component p* 
takes on continuous variations with a distribution. On the other hand, when boosted by 
the Lorentz transformation A, the energy of the ith constituent is given by p® = (Ap*)° in 
a moving composite particle, but their sum, p\ + p®, remains a constant. Even though the 
time-like components of the constituents have these well-defined values which depend on the 
frame of reference, they are often not written out explicitly, for brevity of notation. 

For a proper treatment of spin, we should parallel the treatment in the above sections 
except using Dirac operators instead of Klein-Gordon operators. Alternatively, we can 
adapt the above spinless results to the case of spin by using the fact that we can reduce, 
for two particles, the two-body Dirac equations to Schrodinger-like forms above but with $ 
depending on spin degrees of freedom [Kj, p6| . In a future paper we shall include the spin 
dependent features in more detail. 

In order to deal with multiparticle configurations, we introduce creation and annihilation 
operators of the constituents in a composite particle in its CM frame 



3~T| is then 




(6.23) 



pt,p;) = fet(p*)rft (p * ) | 0) . 



(6.24) 



(We suppress spin, flavor, and color indices.) We assume the general expression 
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{tf(p?)Mpl)} = N(pl)5(p?-pl) 

{<ft(p?), d(p*)} = N(p*)5(p' 2 * - p* 2 ) (6.25) 

and thus 

(p'*,p';\pl,p* 2 ) = N(pl)5(p'* - pl)N(p*)5(p'* - p*). (6.26) 

In the case in which free isolated particles are created and annihilated, one traditionally 
takes either N(p*) = 1 or N(p*) = 2y/p* 2 + m 2 . We emphasize, however, that the above 
momentum in the creation and annihilation operators are not on mass shell but on energy 
shell. In the context of the constraint approach, the individual creation and annihilation 
operators do not produce free-particle states, but rather constituent states within a compos- 
ite associated with a definite total mass and total momentum. Since the aim of this paper 
is a description of the relativistic N— body problem in a separable two-body basis, this is 
plausible. To achieve this, we must determine how the creation and annihilation operators 
will transform under a Lorentz transformation. Let U(A) be our unitary boost operator 
defined so that 

U(A)b\pl)U-\A) = C(p 1A )b\p 1A ) (6.27) 

and 

U(A)d\p* 2 )U-\A) = C(p 2A )rf t (p 2A ) (6.28) 

where p^A is the three- vector part of Ap*. In the case in which free isolated particles 
are produced, the above two conventions lead respectively to either C(p) = \J (Ap)°/p° or 
C(p) = 1. In order to see what these factors become now in the constraint approach, we 
consider 

t/(A){6Vi*), m)}U-\A) = N(pl)5(p'* - pi) 

= C(p 1A )C*(p' 1A )N(p 1A )6(p' 1A - p 1A ). (6.29) 
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To make use of this equation, we need to express the boosted momenta in terms of the 
unboosted momenta. The Lorentz transformation matrix A, defined in Eq. ( |6.11| ), is inde- 
pendent of the momentum of the system being boosted. For the same two-body system, the 
space-like part (Eq. ( |6.13| ) of the Lorentz transformation on two different momenta (for the 
same quark, that is, S\ = e[) yields 

E PP 

PlA - p' 1A = [! + (—- i)_ ](p* - pf). (6.30) 
Using the fact that 5 3 (Ar) = 5 3 (r)/det \A\, we find that 

*(PU - p'ia) = f W - Pi) (6-31) 

so that Eq.( |6.29| ) becomes 

N(pl) = C(pi A )C*(p / 1A )iV(p 1A )^- (6.32) 

We choose the normalization N — 1. Using phase convention with real C, the simplest 
choice is 

This C is associated with the motion of the composite particle. Thus we have 

U(A)¥(pl)U-\A) = /I^(Pia) (6-34) 

and 

U(A)d\ P ;)U- 1 (A) = \/Jrf t (p2A). (6.35) 

(As anticipated above, this contrasts with the on mass shell factor given in standard texts 
(see |32] Eq.(4.2.12)) which refer only to the constituent momenta.) Thus we have 

|M(P) ) = U(A)\M(0)) 

^/^^*,(p;,P^(p; + pa|M(p 1A ),t (P2A )|o). ^ 
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We emphasize that the transformation equations (Eq. (|6.34 ) and Eq. ( |6.35| )) are valid for 



arbitrary p! or p 2 , not just ones that satisfy the rest condition of pi + p 2 = 0. This implies 
that the creation and annihilation operators that go into making up the individual two-body 
interactions will transform in a similar way. (See Eq.( [7.10 ) below.) 



Next we change variables so that the delta function reflects the new total momentum P. 
The inverse of the above Lorentz transformation gives for i = 1, 2 



Pei M .PPp 



Pi = - Z f + P*A + (i " l)^pf^ (6 ' 37) 

which in turn gives 

(Ap , )0 = £ ^ + p_^ (63g) 

Note that with Pia + P2A = P we have 

Pi + P*2 = PlA + P2A - P = 0. (6.39) 

Computing the Jacobian of the above transformation ( 6.37 ) gives 



d 3 p* = ^d 3 p lA . (6.40) 

(Note again how this contrasts with the case of free on-shell particles where d 3 pi = 
d 3 p i Ap i /(Ap i ) .) Hence, 

|M(P)) = —J rf 3 piArf 3 P2A^M(pt,p;)5(P " PlA " P2A)& t (piA)rf t (p 2 A)|0). (6.41) 

The wave function V'mCPijPsO is actually t/>m(Pi(Pia)>P2(P2a)) where the vectors 
{Pi(Pia)> P2(P2a)} are {(pia)a- 1 ) (P2a)a- 1 } • A simple relabeling (not a transformation) 
gives the representation of a general state for a composite particle with momentum P, 

|M(P)> = ^ J d 3 pA<Mp 1A -i, p 2 A-i)5(P - Pi - P2 )6 t (pi)rf t (P2)|0). (6.42) 
In the nonrelativistic limit, this becomes 

im(p)> = / d 3 Pl d 3 p^ M { Vl - 7 p, P2 - m 2 p) 

J [mi + m 2 ) {m 1 + m 2 ) 

x 6(P - Pl - P2 )6 t ( P i)rf t (p2)|0). (6.43) 
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So, each composite state, Eq. ( |6.22| ) with Eq. (|6.42|) , differs from that of the nonrelativistic 



limit by not only replacing Galilean boosts with inverse Lorentz boosts but also by a factor 
of the ratio of the total CM energy to the lab energy. 

In Appendix A we show that the scalar product for two systems in the same internal 
state is given by the three-dimensional momentum delta function times a covariant form 

(M(P')|M(P)) = ^ 3 ( p/ ~ 2 P ) M3 J d*p5(p.P)\i; M (p)\ 2 , (6.44) 

with the wave function having the same dimensions as in the nonrelativistic case. 

Given these preliminaries, we consider how to use this formulation in the calculation 
of meson-meson scattering amplitudes. In this problem one starts with a state | ( 12) (34)) 
consisting of two quark-antiquark states. We model the interaction by the exchange of an 
(effective) gluon corresponding to V. At lowest order, the exchange could not produce a 
final state | (12) (34)) but only | (14) (23)) since the emission of a virtual gluon would leave 
the resultant initial state as two color octet mesons rather than singlet mesons. One would 
thus need to evaluate a typical matrix element of the form 

((14)(23)|y(xi4L)|(12)(34)>. (6.45) 



Inserting / dt\t)(t\ = 1 and using Eq. ( |6.23| ) into the above expression, we can carry out 



the integration in t to obtain a delta function which describes the condition of total energy 
conservation. We have 

f M 12 M l4 M 13 M 2 ^ 2 



({U){23)\V{x li± )\{12)(34)) = 2n6{E 12 + E u -E 13 -E 2i ^ 

x (M(P 14 ), M(P 23 )\V(x ui)|M(P 12 ), M(P 34 )> (6.46) 

where 

|M(P 12 )> = ^ / dVrfW(piA->P2A-)5(Pi2 " Pi - p 2 )& t (pi)rf t ( P2 )|0), (6.47) 



in which P12 is the momentum of the composite with CM energy M\ 2 so that £12 = 
Pf 2 + ^12 an d A^ 1 is the inverse boost to the rest system of the composite. Similar 
expressions appear for |M(P 34 )), |M(P 14 )) , and |M(P 23 )). 
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VII. REACTION MATRIX ELEMENT 



In order to compute matrix elements of the potential, we need its second quantized 
version. We evaluate it in the rest frame of the two interacting constituents and we assume 
that the second quantized form of the potential has the same relation to its first quantized 
form as in the nonrelativistic case. For the reaction A(12) + 5(34) — > C(14) + .0(32), the 
matrix elements for the interaction between a particle and an antiparticle consists of the CI 
and the C2 diagram in Fig. 1. We consider the CI diagram as a representative case. The 
interaction corresponding to this CI diagram takes place between particles 1 and 4 and is 
given by 

V(x u± )\ Pl4=0 = V(± u ) = J d 3 x / 1 d 3 x' 4 d 3 x'ld 3 x'iV(x.' u )5(x.' 1 - x'/)^ - x^lx'^XxX'l 

= J d 3 ^W¥Al<HPi + P4 - pi - p' 4 )v(pi - pD^pi^M^pDHp'i), (7.i) 

in which the integrals include sums over spin, flavor, and color. In this form the indices and 
primes serve to lable the quark color and flavor as well as the momentum. The annihilation 
and creation operators with momentum pi and p[ apply only to particle i. To represent the 
interaction l / (x 14 ), we choose to represent the annihialtion and creation operators in the 
rest frame of the (paticle l)-(particle 4) pair, i.e. the (P14 = 0) frame. In this frame, let 
the total momenta of the two initial meson composite systems be P = P12 +P34, with P i2 
and P34 to be the individual incoming momenta of the two meson composite systems. The 
total energy of the two meson system is given in terms of the respective CM energies of the 
composite particles 

V G = E 12 + O34 = v/p? 2 + Mh + yH + M| 4 . (7.2) 

The total momentum P of the two (composite) particle system is conserved in the scattering 
process so that we can label the matrix element as 

(M(P 14 ), M(P 23 ); P|y(x 14 )|M(P 12 ), M(P 34 ); P). (7.3) 
Let us evaluate the matrix element in the "collider frame" defined by 
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(7.4) 



The Lorentz boost to that frame is given by 



a*-* i ( VWT ~ S n p * Pfc 

Afc - + { ^= i)^ - 

p 

A* — A — - — 

°" l " ^ 



AS^^ (,5, 



and takes us to 



(M(P 14 ), M(P 23 ); P|F(x 14 ) |M(P 12 ), M(P 34 ); P> 
= (M(P; 4 ), M(P C 23 ); 0|f/(A)\/(x 14 )f/- 1 (A)|M(P^ 2 ), M(P§J; 0) (7.6) 



where 



|M(py,M(py ; o 

M 34 



y"rfV^ 3 P2^(PlA r2 i,P2 Al -i)5(Pi 2 - Pi - p 2 )6 t (pi)rf t (p 2 ) 

l rf 3 p3rfW(P3 A - 4l ,P4 A - 4 i)5(P^ 4 " P3 " P4)& t (p 3 )rf t (p4)|0), (7.7) 



(M(P^),M(P§3);0| 

(oi y ^^(Piv^p^o^PM-pi'-p^^piOKp;') 

/ ^V 3 V*(P^' P^)^ - P2 - p£)d(pS)&(p£), (7-8) 



Mm 
M 23 



E23 

and ^ = v/Pf+Mj- 

Note that the integrals in the above state vectors include color summation but not flavor. 
The above matrix element then becomes 

(M(P; 4 ), M(P C 23 ); 0|f/(A)\/(x 14 )f/- 1 (A)|M(P^), M(P§J; 0) 
= (M(P? 4 ), M(Py ; 0| I dV^WVA^Pi + P4 - Pi - p1M Pi - p 4 ) 

x^(A)6t( Pl )dt(p 4 )d(pi)6(p' 1 )^(A)- 1 |M(PS 2 ),M(PS 4 );0>. (7.9) 
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From our earlier arguments, the transformations defined in Eqs. (|6.34 ) and ( 6.35Q are inde 



pendent of the total momentum of the two-body system. Thus, using this, we obtain 

<M(P< 4 ) , M(P C 23 ); 0|f/(A)\/(x 14 )f/- 1 (A)|M(P^ 2 ), M(Py ; 0) 

= (Af(P$ 4 ), M(P C 23 ); 0| J d?pid Z pdh'idh'AVi + P4 - p'i " P 4 )^(Pi " Pa) 

x6t(p 1A )rft (p4A)rf(p ^ )6(p ; A) | M( pc 2);M( pc 4) . 0) ^_ (71Q) 

We can point out that this last factor, due to interaction transformations, is not present in 
the nonrelativistic limit. This matrix element is evaluated in Appendix B, and we find 

(M(Py , M(Py ; 0|f/(A)V(x 14 )f/- 1 (A)|M(P^ 2 ), M(Py ; 0) 
ri/T, M 23 E u M 12 M U 

= -5 3 (P 34 - Pi4 + Pi 2 - p ^)j^Tr-^j^ 

£/ 2 3 M 14 &12 -&34 

X J d 3 Pl d 3 p4lp* D (p 1AA -i, -Pi AA - 1 )V'a((- P 34 + P4a)a 23 1 > -(- P 34 + P4a)a 2 " 1 ) 

Xtfj A ((-P C 23 + P C U - P4a) Ai - 1 ) -(-P23 + P 34 - P4a) A - 1 )?M - P4AA- 1 - P4A A - 1 )^(Pl ~ P*)- i 7 - 11 ) 



From Eqs ( |5.10| ), ( |6.46| ), and ([7.11 ), the reaction matrix element for the interaction Vu is 



hfi[u] 



M 23 M 12 M 34 \ 3 M : 



'14 



J d 3 p 1 d 3 p^ A (p A ,-p A )'ilj B (pB,-pB) 



E 2 3 E\ 2 £"34 / E14 
x ip*c(Pc, ~PcWd(Pd, -Pd)V(Pi ~ Pa), (7.12) 



where 



Pa = (-P^ 3 + P34 - P4a)a-, (7-13) 
PB = -P4AA" 1 ) (7-14) 

Pc = (-Pm + P4a) a ^> ( 7 - 15 ) 

Pd = Piaa- 1 - (7-16) 

Note that the dimension of the wave function leads to an hfi which has the dimension of 
1/mass 2 as required earlier. In the above expression, p^A" 1 ' f° r example, is the space part 
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of (A^Api) and the four energy ratios correspond to the transformations associated with 
the (23), (14), (12), (34) composite particles. The energy ratio for the (14) composite differs 
from the others due to the transformation of the (14) interaction term. The term here 
corresponds to the CI diagram in Fig. I. 

It is interesting to note that compared to the nonrelativistic case, the overlap matrix 
element now involves two major differences. First, the momentum arguments in the wave 
function need to be inversely boosted back to the frame in which the composite particles are 
at rest, as they should be. Second, there are factors of the type M^/E^ appropriate for the 
composite particle in the collider frame. Both effects can lead to substantial modification of 
the magnitude of the reaction cross sections. 

VIII. CONCLUSION AND SUMMARY 

We seek a relativistic formulation of the many-body problem involving both bound states 
and reaction between constituents of composite particles. As a first example, we have 
focused our attention on a system of spinless particles interacting with a scalar and/or 
vector interaction. 

We began by examining the relativistic two-body bound state problem and introduced 
Lagrangian multipliers to write down the most general two-body Hamiltonian. The for- 
mulation using the constraint dynamics allowed a simple separation of the center-of-mass 
and the relative motion. A two-body equation was obtained in the form of a non-relativistic 
Schrodinger equation which connects naturally to the corresponding non-relativistic problem 
in the non-relativistic limit. The two-body equation is independent of the Lagrange multi- 
pliers. The bound state mass is related to the eigenvalue of the non-relativistic problem by 
a simple algebraic equation in the case in which one considers only relativistic kinematics. 
Further relativistic effects show up in this algebraic relationship when the energy dependence 
accompanying the scalar and vector interaction is taken into account. 

For a many-particle system, we considered pair-wise interaction between particles. In 
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constructing the total Hamiltonian, a good choice of the Lagrange multiplier provides a 
simple way to separate the iV-body Hamiltonian into the unperturbed Hamiltonian and 
residual interactions. It also presents a systematic way to use the two-body solution as basis 
states for multi-particle dynamics. The study of the dynamics involves the evaluation of the 
reaction matrix elements of a general two-body interaction in terms of the wave functions 
of the composite particles. 

In rearrangement reactions, because there can be many ways to divide the total Hamil- 
tonian, the evaluation of the reaction matrix elements should not depend on the choice 
of the unperturbed Hamiltonian and basis states. With our formulation, this "post-prior" 
equivalence can be shown explicitly, allowing for a meaningful definition of the perturbation 
expansion and treatment of the reaction dynamics. 

Finally, we give an explicit formula for the reaction matrix elements in terms of the 
composite wave functions. In the relativistic treatment, the important effects include the 
inverse boost of the relative momentum to the frame in which the composite particles are 
at rest, so as obtain the correct wave function. Furthermore, there are factors of M^/Ey 
in the collider frame for the composite particles. These relations will be useful when we 
apply the present formulation to many problems in nuclear and particle physics such as 
meson-meson scattering. An aim would be to see how this approach modifies the results of 



the nonrelativistic formalism as present in [18 



The results we have obtained are very encouraging. We should in future work carry out 
a calculation for the relativistic I = 2 tttt scattering, to compare with the non-relativistic 



results of Barnes and Swanson [18[ and with experimental data. We should also extend our 
considerations to include the spin degree of freedom and more complicated interactions in 
the constraint description. 
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APPENDIX A: SCALAR PRODUCTS 



Let us consider the general scalar product (M(P')|M(P)) 

( M(P')|M(P)> = ~J rfVlA^V 2 A'(OM(p , 2A0KPlA0V'M(p , l,P2)5 3 (P / -p / lA' -P2A') 
X J d 3 p 1A d 3 p 2 A4>M(Pl,P2)5(P - PlA - P2A)& t (PlA)^ t (P2A)|0) 

= ^ j d 3 p' 1Al d% A , J d 3 PlA d 3 p 2A r M (p'l, PsMPi, P2)5 3 (P' - P)5 3 (P - PlA - P2A) 

x5 3 (piA-p' 1A ,)5 3 (p 2A -P2A'), (Al) 

where E = VP 2 + M 2 , £" = VP' 2 + M 2 . 

Now the total momentum delta function makes A = A', which in turn implies that 

the two delta functions that come from the creation and annihilation operator force the 

arguments of the two wave functions to be the same. Thus 

M 2 r 

( M(P')|M(P)> = 5 3 (P' - P)_ J rfW 3 P2A^M(Pl,P2)^M(Pl,P 2 )5 3 (P " PlA " P2A) 

M 2 r 

= 5 3 (P' - P)— J rf 3 pirfWM(PlA-i,P2A-i)^M(plA-i,P2A-0^(P ~ Pi ~ P2). (A2) 

On the right-hand side let 

P = Pi + P 2 , 

P = M Pl " M P2 ' (A3) 

so that the integral becomes 

/ d*Pd 3 P r M ((^P + p)a-i, ((§ P - p))a- 1 )^((§P + p)a-i, ((§ P - p))a-0*(P " P) 
= / ^m((§P + p)a-i, ((ff P - p))a- 1 )Vm((§P + p)a-i, ((ff P - p'))a-0 
= / ^m((0 + p')a-i, ((0 - p))a-0^m((0 + p) A -i, ((0 - p)) A -i) 

= J d 3 p|VM(PA-0| 2 = / rfV(p°)I^M(A-V)| 2 . (A4) 

But p° = p ■ A _1 P/M. So taking p = Ap and using d 4 p' = rf 4 p, we have the following 
manifestly covariant scalar product 

/ d 4 pS(p°M M (A^p)\ 2 = M J d 4 p5{p-P')\4> M {p)\ 2 . (A5) 
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APPENDIX B: EVALUATION OF THE REACTION MATRIX ELEMENT 

In this Appendix we evaluate the matrix element of Eq. (|7.10|) , 



(M(Py , M(P C 23 ); 0|f/(A)K(x 14 )f/- 1 (A)|M(P^ 2 ), M(P U ); 0) 
= (M(Py , M(P C 23 ); 0| J d 3 p';>d 3 pl>d 3 p[d%5(p';' + p'i' - p[ - p 4 )V( Pi " - p 4 ) 

xf/(A)6t(p' 1 ")rf t (pDrf(pl)Kp , i)^(A)- 1 |M(P? 2 ),M(P^ 4 );0) 
= (M(Py , M(P C 23 ); 0| J dhWvW.dhW" + pi' ~ P'l " P',MP7 ~ P 4 ) 



in which the integrals include sums over flavor and color. This requires us to compute 



(0|d( P :)Kp;Orf(pDKP^&^P;A)rf t (P4A)rf(P4A)Kp / lA)& t (Pl)rf t (P 2 )& t (P3)rf t (P4)|0) 

= (0|[*(p& - pD5(Pia - PD^PDKPi) - ^(P4A - P4>(Pia - PD^PsOKPx) 

-^(P4A - PSWPIA - Pl)d(P4)6(Pa) + *(P& - Pl'WPlA - PlMP3>(P 2 )] 
X [<KP 4 A - P2)*(PlA - Pl) fot (P3)^(P4) - $(P4A ~ P^PlA ~ Pl^M^M 

-S(p' 4A - p 2 )6(p' 1A - p 3 )k t ( P iM t (P4) + 6(p' 4A - p 4 )5(p' 1A " P3)& t (pi)rf t (p 2 )]|0) (B2) 



in which the delta functions include flavor and color indices. If we assume that flavors for 1 
and 2 are distinct from those of 3 and 4, then of the sixteen terms above the only one that 
survives is 



(o|*(p& - pMp? a - p1)S(p' 4A - P4 )<Kp' 1A - pi)*(pS - P3 )<H P2 - P2 )|o). (B3) 



<M(P? 4 ) , M(Py ; 0|f/(A)nxi4)^ 1 (A)|M(Py , M(P C U ); 0) 
= - / WrfWplW^^ - p'x'a - p2a)^(pV-i, P3A-0 

j/ 14 14 2o 

X<K P 2 3 - P 2 - P'DMp't^-i, P2 A -l)^(P?2 - P 2 - Pia)V'b(P3A- 4 1 ) P4AA 3 - 4 1 )^( P 34 ~ P 4 A ~ Ps) 



&t(p-)dt (p - )d(p ^ )6( p' iA) | M( pC 2))M( pC 4) . ) 




(Bl) 



Thus (using the notation = 



P?? + M?-), we obtain 



+p^-pi- P i)v( P r-Pi) 
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(B4) 



-E 2 3 Afi4 -E12 -E": 



'34 
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We perform four of the remaining six volume integrals of the first set of integrals. In 
particular, we perform integrations over (Pp'^d^p'^p'^p 1 ^- We use the first delta function 
p'l' A = — p'" A , the second gives p' 2 ' = P 2 3 — P3, the third gives p' 1A = ~P\ 2 — p 2 and the 
fourth gives p 3 ' = P34 — p 4A . The argument of the remaining delta function is then 

Pi' + P4" ~~ Pi — P4 = Pua- 1 ~ P12A- 1 + P34A- 1 ~~ P23A- 1 = P34 — Pi4 + P12 — P23) (B5) 

which corresponds to overall momentum conservation in the frame in which we evaluate the 
matrix element (the P14 = frame). Then, the reaction matrix element becomes 

(M(P< 4 ) , M(Py ; 0\U(A)V(x l4 )U-\A)\M(P c l2 ), M(P C U ); 0) 
= -5 3 (P 34 - Pi4 + P12 - P23) / dWri^P^-i, ( p i4 - p'ia)a-) 

X^(P C 23 - P34 + P4a)a-, (P34 - P4a)a 2 -0^((P? 2 " P° 23 + P34 - P4A)A-, (P23 " ?34 + P 4 a)a"0 
/ //t>c / \ / \t>/ "I I \ ^23 ^14 M12 M34 

x^((P 34 - P4a)a- ? P4Aa 3 -0^(Pi - p 4 )_ — — — . ( B6 ) 
But Pi4 A -i = = P12A- 1 = P34A- 1 = P 23a-i so that 
(M(P< 4 ) , M(P C 23 ); 0|f/(A)\/(x 14 )f/- 1 (A)|M(P^ 2 ), M(P< 4 ); 0) 

= -5 3 (P 34 " Pl4 + Pl2 - P 23 ) / ^'/'^^(p'/aa-, " P'lAA- 4 0#((-P34 + P^A^ (P34 " PLOa'O 

XV(P'1' - Pi)lM(-P§3 + P 34 - P4a)a-, (P23 - P34 + P4a)a-)^( " P^i, P^i ) ^ ^ , 

23 14 12 34 

which corresponds to the indicated amplitude in which a gluon is exchanged between particle 
1 and particle 4. This produces the form in the text (where we dropped the primes on the 
two integration variables). 
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